Checkerboard superconducting order and 
antinodal Bogoliubov quasiparticle interference 



V. I. Belyavsky, V. V. Kapaev, and Yu. V. Kopaev 
P. N. Lebedev Physical Institute of Russian Academy of Sciences, Moscow, 119991, Russia 

Numerical study of momentum-dependent gap function is presented to make clear the origin 
of superconductivity in copper oxides. We claim that antinodal region with pronounced nesting 
feature of the Fermi contour gives rise to superconducting pairing with large momentum under 
screened Coulomb repulsion. Such a pairing results in both spatial checkerboard pattern of the 
superconducting state below T c and a gapped state of incoherent pairs in a broad temperature range 
above T c . We explain the momentum dependence of the coherent spectral weight detected in angle- 
resolved photoemission spectroscopy and predict antinodal Bogoliubov quasiparticle interference 
other than observed in the nodal region. 
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I. INTRODUCTION 

Angle resolved photoemission spectroscopy (ARPES) 
of underdoped cuprates at temperatures exceeding con- 
ditional upper boundary T* of the pseudogap (PG) state 
evidences in favour to large simply connected Fermi con- 
tour (FC) typical of the conventional Fermi liquid. How- 
ever, below T* , the FC is seen as transformed into dis- 
connected arcs disposed in the nodal regions. Cooling 
from T* down to superconducting (SC) transition tem- 
perature T c results in a decrease of arc length down to 
zero. Thus, the FC degenerates into four points that give 
rise to the nodes of d - wave SC order parameter arising 
below T c . It seems fairly natural to conclude that, within 
the framework of d - wave pairing concept, the SC order 
parameter has its maximal value exactly in the antinodal 
directions, 1 

For this reason, it might seem quite probable that 
low-temperature properties of d - wave superconductor 
should be determined by low-energy quasiparticle exci- 
tations only in the nodal region of the momentum space 
that is in vicinities of the SC gap nodes on the diag- 
onals of the Brillouin zone. Taking into account that 
Bogoliubov quasipatricle interference (QPI), observed in 
the nodal region, disappears near the end points of the 
FC arcs^ one might lead to a conclusion that only the 
nodal region gives rise to superconductivity whereas the 
gap observed in the antinodal region should be attributed 
to an incoherent PG stated However, a coherence in the 
antinodal region becomes apparent both in the ARPES 
study 5 and also in the Andreev St James experiments. 4 
Therefore, in spite of the fact that high-energy QPI is not 
detected for the present, one can believe that the antin- 
odal region should contribute a coherent state as well. 

We have argued^ that both PG and SC states arise 
exactly in the antinodal region with pronounced nest- 
ing of the FC as spatially inhomogeneous incoherent and 
coherent states of pairs with large momentum, respec- 
tively. The nodal region gives rise to conventional SC 
pairing with zero momentum which, together with the 
pairing with large momentum (/^-pairing) in the antin- 



odal region, forms a biordered SC state in the whole of 
the Brillouin zone. 

Kinematic constraint, inherent in it-pairing in the 
antinodal region, can result in oscillating real-space pair- 
ing interaction. Indeed, momenta of both particles com- 
posing SC pair with nonzero total momentum K should 
be either inside or outside of the FC. For this reason, a 
set of one-particle states turns out to be kinematically 
excluded because of the fact that such states cannot con- 
tribute into the states of it-pairs. It means that any 
scattering between such excluded states should be for- 
bidden when one defines the interaction leading to a rise 
of a bound state of it-pair. An exclusion of a set of the 
Fourier components from the screened Coulomb inter- 
action results in the fact that corresponding real-space 
it-pairing interaction exhibits an oscillation outside of 
small-distance repulsive core as shown schematically in 
Fig. 1. It should be noted that there is an analogy be- 
tween this oscillation and well-known Friedel oscillation 
that arises owing to Kohn singularity of screening en- 
hanced by nesting of the FC. 

Besides the fact that two-particle problem with os- 
cillating potential leads to a bound state of the rela- 
tive motion of it-pair j£ it can also produce a quasi- 
stationary state (QSS}£ similar to the Gamov's state 
of alpha-radioactive nucleus^ SC gap function A sc (fc), 
depending on relative-motion momentum k of it-pair, 
as a solution to the mean-field self-consistency equation, 
arises as a result of the instability of the ground state of 
the normal Fermi liquid with respect to a rise of it-pairs 
in the bound state. This function can be expressed in 
terms of Gorkov's anomalous averages describing SC con- 
densate of it-pairs. Due to a phase coherence of the SC 
ground state, these averages become nonzero below T c . 
QSS with positive energy, following from the two-particle 
problem, can be considered as an evidence in favor to one 
more instability of the Fermi liquid. We associate such 
an instability with incoherent it-pairs existing above T c . 
It means that Gorkov's anomalous averages and corre- 
sponding gap function A pg (k) vanish under averaging 
over phases of it-pairs but mean square gap function 
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FIG. 1: Real-space pairing potential U(r) (schematically). 
Energies _E 4 and E q correspond to bound and quasi-stationary 
states, respectively. Barrier height E b corresponds to a break 
of the pair without tunnelling through the barrier. 



remains nonzero up to T* according to the hypothesis 
advanced by Emery and Kivelson^ 

Thus, one can conclude that, in the temperature range 
from T* down to T c , iiT-pairs form incoherent PG state 
as off-diagonal short-range order (ODSRO). Off-diagonal 
long-range order (ODLRO) arises as SC condensate of co- 
herent K-pairs below T c . Both ODLRO and ODSRO 
states can be described in terms of Gorkov's Green 
functions^ A phenomenological BCS-like form of the 
coherent contribution to the normal Gorkov's function 
can be written as 



G(w;k) = z(k) 
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where E(k) and 2u±(k) — 1 ±£(k)/E(k) are quasiparti- 
cle energy and coherence factors, respectively, 



2£(fc) = e(K/2 + k)+ e(K/2 - k) 



(1) 



is the kinetic energy of the K - pair of particles with mo- 
menta K/2 ± k, e(k) is electron dispersion with respect 
to chemical potential /i and z{k) is quasiparticle weight. 
Two terms in G(uj;k) can be referred to .K"-pairs above 
and below the FC, respectively. Diagonal Green func- 
tion G(uj;k) describes ODSRO state corresponding to 
the existence of non-coherent QSS of K - pairs above T c . 
Transition from the bound paired state into long-living 
QSS corresponds to small but finite decay T — T(u>; k) 
whereas transitions into stationary states above barrier 
energy E b (Fig. 1) should be associated with an infinites- 
imal decay, 7 — > +0, leading to conventional Fermi-liquid 
behavior of G(uj;k) above T*. Thus, a rise of QSS re- 
sults in a non- Fermi-liquid behavior of G(u>; k)) that can 
be related to the PG state. 

The SC state below T c should be described by both 
normal and anomalous Gorkov's functions. Taking into 
account the fact that PG function A (fe), averaged over 
random phases, vanishes whereas A sc (fe) ^ below T c , 
one can introduce anomalous Gorkov's function F + (u>; k) 



F+(uj;k) = -z(k) 
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Such an approach directly leads to uniform descrip- 
tion of both SC and PG states in underdoped cuprates. 
One can see that repulsive Coulomb pairing in the antin- 
odal region necessarily results in rather complicated mo- 
mentum dependence of the SC gap and PG functions, 
A sc (fc) and A pg (k), with energy scale e ~ I eV of their 
domains of definition 11 in contrast to considerably less 
scale of about Debye energy e D that arises in the case of 
phonon-mediated SC pairing. It is very likely that the 
high energy problem^ arising, in particular, in the op- 
tical conductivity of the cuprates ^ might be associated 
with high energy scale of the antinodal K -pairing. We 
believe that _RT-pairs are the main players in the high- 
temperature superconductivity of the cuprates. 

Recently, Tsvelik and Chubukov 14 considered SC pair- 
ing on mutually orthogonal pairs of perfectly nested seg- 
ments of the FC in semiphenomenological way. They 
presuppose that one-dimensional SC order arises only on 
these segments coupled with a momentum-space Joseph- 
son links to give rise to two-dimensional superconductiv- 
ity. Actually, it is implicitly supposed that such SC state 
can arise owing to SC pairing with nonzero momentum. 
Also, it is supposed that SC order with the same mo- 
mentum is induced on the rest unnested part of the FC 
by the order on the nested segments due to a proximity 
effect in the momentum spaced so that this induced order 
cannot penetrate deep into the nodal region. It should 
be emphasize that such a model^ differs essentially from 
the biordered SC stated 

In this paper, we study the mean-field If -pairing prob- 
lem numerically to fall outside the weak coupling lim- 
its employed in our previous approach to the K -pairing 
problem We show that the SC gap function with a 
nontrivial nodal line corresponds to a checkerboard pair 
density wave (PDW) SC state and results in fairly natu- 
ral explanation of the angle dependence of a partial sup- 
pression of the coherent spectral weight in the antinodal 
region observed by Kondo et al£ We believe that QPI, 
other than observed in the nodal region,— could be de- 
tected in the antinodal one as well. We also show that Im- 
pairing can originate spatial checkerboard pattern with- 
out any driving insulating order in contrast to a scenario 
of a rise of PDW coexisting with a charge density wave 
(CDW)^ 



II. K - PAIRING PROBLEM 



In the case of JC-pairing, the gap function is defined 



A(fe) = ^t/(fe,fe')(< 
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where U(k, k') is screened Coulomb interaction matrix 
element, operator c ka annihilates electron with momen- 
tum k and spin polarization a. Anomalous average in 
Eq. @, describing SC condensate of K - pairs, becomes 
nonzero below T c . The gap function should be a nontriv- 
ial solution to the self-consistency equation, 
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Here, n(k) = (e E ( k ^ T + 1) 1 is a quasiparticle occupa- 
tion number and quasiparticle energy has the form 



E(k) = n(k) ± y/£ 2 (k) + A 2 (k) 
2r?(fc) = e{K/2 + fe) - s(K/2 - k). 



(4) 
(5) 



It should be noted that, since e(— k) = e(k) owing to the 
time-reversal symmetry of the dispersion relation, quasi- 
particle spectrum ^ turns out to be gapped on the whole 
of the FC in the case of pairing with zero total momen- 
tum. In the case of _KT-pairing, it can be gapped only 
on those parts of the FC where \r](k)\ proves to be small 
enough. In addition, the quasiparticle spectrum becomes 
asymmetrical with respect to /i. 

Summation in Eqs. and §5§ should be performed 
over all momenta of the relative motion which can form 
pairs with given total momentum K. One can see that 
these momenta belong to a K - dependent domain of 
the momentum space (domain of kinematic constraint) 
because of the fact that the momenta of both particles 
composing SC pair with given total momentum should be 
situated either inside or outside of the FC. This means 
that some part of the momentum space turns out to be 
excluded from the sums in Eqs. (|2I3[) . 

Since the kinetic energies of both particles composing 
SC pair with K = can be equal to /i, the low-energy 
limit in the sum §5§ corresponds to £ = whereas the up- 
per limit is formally restricted by a half-width of the con- 
duction band of the order of \i. However, in the Bardeen- 
Cooper-Schrieffer (BCS) theory^ such upper limit (De- 
bye phonon energy e D ) appears as an energy scale of 
a layer enveloping the FC where electron-electron scat- 
tering results in an effective attraction between electrons. 
As a result, pairing interaction energy U(k, k') in Eq. ([3]) 
can be qualitatively associated with an effective coupling 
constant V* that can be estimated as^I 



V* = V - 
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where g is density of states per spin, V is a pairing con- 
stant due to electron-phonon interaction defined inside 
the layer, U is average Coulomb energy. Thus, in the ef- 
fective pairing constant, Coulomb repulsion appears with 
a logarithmic weakening. In the case when V*g <C 1, the 
mean-field approach results in a conventional BCS energy 
gap, 
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FIG. 2: Near nested FC (solid line) corresponding to electron 
dispersion Eq. © typical of the cuprates. Dashed lines rep- 
resent isolines close to the FC, numbers near the isolines are 
electron energies according to Eq. (0 with t = 0.5, t' — —0.15, 
t" — 0.07 eV. Here, K is total momentum of K-pair, Q is 
nesting momentum. Length L of near rectilinear segment of 
the FC is shown at given mean square energy deviation S (the 
width of the shadowed strip) of the FC from the rectilinearity. 



that appears in consequence of the logarithmic singular- 
ity of the right-hand side of Eq. ([3]) . This singularity is 
primarily formed in an energy range near the low-energy 
limit, therefore, extension of this range might lead to a 
progressive accumulation of the singularity along with 
the formation of a non-singular (regular) contribution 
into Eq. Q . One can treat the preexponential in Eq. ([7|) 
as a characteristic energy scale beyond which the non- 
trivial solution to the self-consistency equation becomes 
weakly sensitive to the upper limit. 

All these speculations can be referred to the K - pair- 
ing problem. However, in such a case, the logarithmic 
singularity becomes apparent if and only if kinetic en- 
ergy of K - pair vanishes not at isolated points, as it 
were most likely in the case of arbitrary FC, but on finite 
pieces of the FC on which mirror nesting condition, 



s(K/2 + fe) - e{K/2 - k) = 0, 



(8) 



A = 2 £D exp(-l/y* 5 ), 



(7) 



should be fulfilled at given K . One can see that, for a 
suitable K, this condition can be fulfilled in the case of 
rectilinear parallel segments on the opposite sides of the 
FC. It is obvious that K should be directed along these 
segments. 

FC, typical of the cuprates, and isolines close to it can 
be described satisfactorily by electron dispersion 

e(k x ,k y ) = t — 2t(cosk x + cosk y ) — At' cos k x cos k y — 
- 2t"(cos2k x + cos2k y ) (9) 

with fitting parameters t Q = 2eV t — 0.5 eV, t'/t = 
—0.3, t" = 0.14. Here, k x and k y (in units of 7r/a; a is 



interatomic distance) are momentum components corre- 
sponding to the antinodal directions. 

One can choose one of the coordinate axes (fc ) along 
K directed parallel to antinodal near rectilinear segments 
of the FC as shown in Fig. 2. Then, with preassigned ac- 
curacy 5, £(fe) < 5 if k y corresponds to near rectilinear 
segment of the FC. Therefore, the singular contribution 
into turns out to be proportional to the length L of 
such a segment. Summation over the other component 
(k x ) leads to an accumulation of the singularity, however, 
in contrast to the case K — 0, a gradual deviation from 
the FC results in a progressive increase of the difference 
between the kinetic energies of the particles composing 
K - pair. This leads to increasing deviation from mirror 
nesting condition (|8|), so that, finally, the accumulation 
turns out to be completed when k x attains a value corre- 
sponding to energy scale e much lesser than fi. It should 
be noted that, in the case of K - pairing, e appears as 
generic energy scale originating from mirror nesting fea- 
ture of electron dispersion. This scale should be related 
to a preexponential of the gap function in the case of 
small effective coupling constant. Thus, one can con- 
clude that the nontrivial solution to the self-consistency 
equation should be weakly sensitive to the part of the 
momentum space corresponding to e > e Q . 

Since equality ([H]) is fulfilled only approximately, one 
can choose length L in order that mean square deviation 
of the FC from the rectilinearity were less than a pre- 
assigned value corresponding to energy scale S. Strictly 
speaking, nonzero 6 eliminates the singularity because of 
a rise of a lower limit cutoff in the sum Eq. ([3]). Similar 
cutoff appears in the Fulde - Ferrell - Larkin - Ovchin- 
nikov (FFLO) problem of SC pairing with small total pair 
moment um i 1 8 ' 1 9 Therefore, nontrivial solution to equa- 
tion can exist if the effective coupling constant ex- 
ceeds certain S - dependent value. Magnitude A' of the 
gap function can be roughly estimated as 



A' = VA(A - 26), 



(10) 



where A is the magnitude corresponding to perfect mir- 
ror nesting that is to exactly rectilinear segment of the FC 
with length L. Positive function A'(<5) has a maximum 
at certain S = S m < A/2. Indeed, A' — > if 5 — > (then, 
generally speaking, L —> 0, so that there is no singularity 
in the self-consistency equation: A — > 0). At 26 > A, 
the magnitude of the gap function vanishes, therefore, 
a maximum value of A' exists at < 6 m < A/2. A 
choice of length L of near rectilinear segment of the FC 
at given 5 predetermines total pair momentum K. It 
is clear that, because of kinematic constraint, the abso- 
lute value of K coincides with L/2, as one can see from 
Fig. 2. Since maximum value of A' corresponds to 6 m , 
the absolute value of the momentum of K - pairs in the 
SC condensate should be taken as K — L(S m )/2. Vari- 
ation of the FC with doping x in hole doped cuprates 20 
results in a conclusion that K should be dependent on 
x. Note that there is no contradiction between such a 
dependence, following from dispersion @, and doping 
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FIG. 3: Calculated momentum dependence of the antinodal 
gap function (gradation of grey, in meV) shown in a region of 
the Brillouin zone that primarily forms the singularity of the 
self-consistency equation. White curves represent the nodal 
line structure of the gap function, dashed line is the FC corre- 
sponding to that shown in Fig. 2. White regions, in which the 
nontrivial solution is absent due to the kinematic constraint, 
appear under shifting of the FC at ±K/2 along k y (inside 
these regions, there are no one-particle states that could form 
a pair with total momentum K — 0.2 -k j a) corresponding to 
a maximum magnitude of the gap function. 



dependence of spatial periodicity of checkerboard PDW 
seen in tunnel data.— 

Comparatively small vicinity with energy scale e of 
the strip with length L(5 m )/2 and width corresponding 
to energy scale 5 m can be considered as the region of 
the momentum space that primarily forms the singular- 
ity of the self-consistency equation. Following Ref. [10], 
one can renormalize the kernel of this equation and re- 
duce Eq. ([3]) to a sum over momenta belonging to such a 
vicinity only. Renormalizcd kernel, defined in this vicin- 
ity, can be written as^ 2 . 



W{k,k') 



K + .9 ln OAo)' 



(11) 



and can be treated as a pairing pseudopotential cor- 
responding to oscillating real-space pairing interaction. 
Here, </>„(&) and A„ are eigenfunctions and eigenvalues of 
kernel U(k, k'), respectively— We believe that the vicini- 
ties with energy scale e of the antinodal near rectilinear 
segments of the FC include electron states that mainly 
contribute to scattering resulting in _RT-pairing. 
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Since characteristic sizes of the vicinity are much less 
than characteristic Fermi momentum, region of attrac- 
tion in the real space proves to be more deep and ex- 
tended with respect to that due to Friedel oscillation. 
Such oscillating interaction can provide both bound state 
and QSS of the relative motion of K - pair. In the 
mean-field approach, the bound state appears in tem- 
perature range < T < T c as nonzero anomalous aver- 
ages, (c K /2-ki^K ^ ^' *^ a ^ determine gap func- 
tion Eq. (J2]). It should be noted that, in the case of small 
K (for example, in the FFLO state), real-space oscil- 
lation of the pairing interaction becomes weak enough 
because of considerable extension of the corresponding 
vicinity forming the singular contribution into the self- 
consistency equation. 

To study iiT-pairing problem numerically, we use a 
step-wise approximation of the pairing interaction^ as- 
suming that pseudopotential (fITj) has a constant value of 
about 10 eV inside a vicinity of near rectilinear segments 
of the FC. Energy scale of such a vicinity is determined 
from the above mentioned condition that calculated gap 
function magnitude should become actually independent 
of this scale beginning with certain s . 

Numerical study of Eq. |3]) at T = reveals highly 
complicated momentum dependence of gap function 
A(fc), shown in Fig. 3, with a few closed nodal lines 
crossing the FC. Topological feature of the gap func- 
tion, shown only inside the part of the Brillouin zone 
that primarily contributes into the singularity of the self- 
consistency equation, turns out to be weakly dependent 
on small variation of the parameters of electron disper- 
sion and magnitude of pairing interaction. According to 
rough estimation following from (|10[) , maximum value of 
K - dependent magnitude of the gap function can be as- 
sociated with K close to 0.2 7r/a. Domain of definition of 
the pairing pseudopotential includes all energies e < e , 
where e is relative to a distance between the FC and 
the boundary of this domain. As follows from numerical 
solution to Eq. a gradual decrease of the momentum 
corresponding to upper limit k r in the sum over k x in 
the self-consistency equation with renormalized kernel, 
at first, does not affect the magnitude of the gap func- 
tion. Then, beginning with certain value of k r , that can 
be associated with a boundary of the domain of definition 
of W(k, k'), the magnitude tends to zero with a decrease 
of k r . This gives a possibility to determine energy scale 
e ~ 0.3 eV of this domain forming the singularity of the 
self-consistency equation. 



III. CHECKERBOARD ODLRO 

Visualization of a checkerboard PDW— can be consid- 
ered as an indirect evidence in favour of the fact that such 
a state originates from nesting feature of the FC typical 
of the cuprates. It should be noted that near rectilinear 
segments on the opposite sides of such FC ensure not only 



mirror nesting condition ([8]) but also nesting condition 

e(k + Q) + e(k) = (12) 

at certain nesting momentum Q which, in general, is in- 
commensurate with total momentum K of SC pair as 
shown in Fig. 2. Under condition (fT2"| . the logarithmic 
singularity can arise in an insulating pairing channel that 
gives rise, for example, to CDW. Such an insulating pair- 
ing can compete or coexist with K - pairing in a way 
considered a long time ago in the case of the coexistence 
of conventional (K = 0) SC state and CDW^i 

One can compare efficiencies of both channels with the 
help of a crude estimation of lengths L K and Lq of near 
rectilinear segments forming singularities in the SC and 
insulating channels, respectively. These lengths, at given 
5, the same in both channels, can be found from inequal- 
ities 

\e(K/2 + k)-e(K/2-k)\ < S, 

\e(k + Q) + e(fe)| < 5, (13) 

selecting the regions in the momentum space in which 
mirror nesting or nesting condition, respectively, is satis- 
fied with preassigned accuracy. If boundaries of these re- 
gions intersect the FC, lengths L K and Lq should be de- 
fined as distances between the corresponding intersection 
points. Both pair momentum K and nesting momentum 
Q should be selected in a way to ensure maximum val- 
ues of corresponding lengths L K and Lq, respectively. 
Momenta K and Q depend on a form of the FC vary- 
ing with doping. Therefore, interrelation between L K 
and Lq varies with doping as well. Calculated variations 
of L K and Lq with doping are shown schematically in 
Fig. 4. A comparison of L and Lq shows that, in the case 
of electron dispersion Eq. ([5]), nesting dominates mirror 
nesting in electron doped compounds. On the contrary, 
the opposite case of hole doping gives an opportunity of a 
rise of such a range of x where mirror nesting dominates 
nesting. 

In this range, SC order arises due to K - pairing 
whereas insulating order (CDW with Q = 2K) with low 
spectral weight can exist as induced by the correspond- 
ing PDW and, therefore, turns out to be hardly detected. 
Conversely, dominating insulating order in the form of 
CDW can induce the PDW— as a modulated superfluid 
density which, due to low spectral weight, seems to be un- 
detectable. Thus, one can conclude that, if mirror nest- 
ing dominates nesting, SC state in the form of PDW of 
K - pairs can exist without any driving insulating order. 

It should be emphasized that nesting feature of the FC 
can lead to a giant enhancement of the singularity in both 
pairing channels i 24 i 25 In the insulating channel, nesting 
of the FC transforms weak Kohn anomaly into the CDW. 
In the SC channel, mirror nesting gives rise to K -pairing 
which, owing to kinematic constraint, ensures extended 
and deep oscillation of real-space screened Coulomb pair- 
ing interaction and, consequently, results in T c consid- 
erably greater than following from Kohn-Luttinger SC 
pairing 2 ^ with K = and angular momentum 1^0. 
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FIG. 4: Competition between nesting and mirror nesting: 
doping dependence of lengths Lq and L of the optimal near 
rectilinear segments of the FC under nesting and mirror nest- 
ing conditions, respectively. 



Thus, -pairing leads to an independent order so that, 
in such a case, there is no need to take into account a co- 
existence of SC and insulating ordered states^ to invoke 
spatially inhomogeneous SC state in the form of PDW. 

One can define anomalous averages {ip ^{r')il>^(r)) in 
the real space that corresponds to the momentum-space 
anomalous averages, (c K / 2 -k{^K /2+ky '■> ar i sm g due to 
.Ff -pairing. Here, fermion field operator tp a (r) annihi- 
lates electron with spin polarization a and radius vector 
r. Nonzero anomalous averages can be considered as an 
order parameter corresponding to _K"-pairing. In the case 
of the two-dimensional C 4 orbital symmetry, there are 
four crystal equivalent pair momenta Jf •, j = 1,2,3,4. 
Therefore, real-space representation of the order param- 
eter should be written as a superposition 

4 

(^(r')^ T (r)> V 3 {P) ■ expQKjR), (14) 

3=1 

where R= (r + r')/2, p = r — r' are center-of-mass and 
relative motion radius vectors of Jf-pair, 



(15) 



can be considered as a real-space wave function of the 
relative motion of JC-pair. Here, N is a number of unit 
cells of the two-dimensional system, summation over k 
should be performed inside the domain of kinematic con- 
straint corresponding to each momentum Kj. Coeffi- 
cients jj, corresponding to SC state, should be deter- 
mined by one of the irreducible representations of the 
symmetry group C 4 . A choice of the irreducible repre- 
sentation establishes the orbital symmetry of the order 



parameter. Since 7 X = — 7 2 = 7 3 = —74 in the case of 
rf-wave orbital symmetry, a checkerboard spatial pattern 
of the order parameter follows from Eq. (fT4"|) immedi- 
ately. One can see that <i-wave order parameter (fLt]) 
corresponds to a current-less SC state, therefore, in this 
respect, it is similar to Larkin - Ovchinnikov immobile 
wave solution^ in contrast to Fulde - Ferrell running 
wave^ of the FFLO problem. 

Nonzero anomalous average (c K / 2 -ki^K /2+k]) a P~ 
pears as a result of averaging of the product of two an- 
nihilation fermion operators over the canonical ensemble 
in which total particle number N fluctuate with respect 
to certain mean value In such an ensemble, all of 

the states with different N close to N should be coher- 
ent so that pair correlation function Eq. (fT4)) describes 
ODLRO of K - pairs in the SC condensate. Above T c , 
phase coherence of the ground state turns out to be lost 
due to the fact that iiT-pairs in the states with differ- 
ent N have got random center-of-mass phases. Owing to 
above-mentioned instability of the ground state with re- 
spect to a rise of QSS of if-pair, relative motion phase of 
the wave function of such a pair, included into coefficients 
7,-, can remain locked up to temperatures far above T c . 

One can think that a loss of relative-motion phase co- 
herence with heating might go through two steps. At 
first, d-wave current-less superposition (| 14[) can be de- 
composed into two orthogonal dimer superpositions with 
7i = ±73: 7 2 = 74 = and 7 2 = ±j 4 , 7l = 7a = 0. After 
that, at greater temperature, dimer state can be disinte- 
grated into free _KT-pairs which survives up to their break 
at a temperature that can be associated with the upper 
boundary of the PG state. Temperature range, corre- 
sponding to lost center-of-mass phase coherence but sur- 
vived relative-motion phase coherence, can be referred 
to the region of the PG state in which off-condensate SC 
pairs can appear as spatially inhomogeneous ODSRO. 
We believe that spatial pattern, like that observed by 
Kohsaka et al.,- is described by current-less superposi- 
tions Eq. (|14p in which coefficients 7 ■ correspond to ran- 
dom dimer configurations. 

Recently, Berg et ali^l have considered dimcr-like 
( "striped" ) ODLRO in the framework of the concept of 
SC pairing with large momentum. Note that, as follows 
from numerical study of Hubbard model on 4 x 4 square 
lattice^ d-wave checkerboard order as the ground state 
seems to be favorable with respect to dimer-like one. 

Momentum dependence of ODLRO parameter A sc (fc) 
determines the angle dependence of the spectral weight, 
W CP ((j)), of the SC coherent peak appearing in the 
ARPES spectra below T c . In the case of d-wave super- 
conductor, it is tacitly assumed that SC order parameter, 
taken on the FC, is proportional to cos20, where Fermi 
angle <fi is polar angle in the momentum space counted 
from the antinodal direction. Therefore, W CP ((/>) should 
be a monotone function in the angle range < (f> < it/ A 
between the antinodal and nodal directions. However, 
the ARPES study^ shows unambiguously that the SC 
spectral weight turns out to be highly non-monotonic: 
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at first, W CP (4>) increases from zero at <j> — 7r/4 up to 
a maximum at certain <j) m and then exhibits a consid- 
erable decrease if cj> — > 0. The spectral weight in the 
antinodal region is strongly dependent on doping. Such 
a non-monotonic behavior of W GP ((f>) is explained 3 by a 
competition between superconductivity and an insulating 
state developing in the antinodal region with pronounced 
nesting feature of the FC. The insulating state should 
result in a depletion of the SC pairing channel and, in 
consequence of a decrease of the SC order parameter, 
in a lowering of the spectral weight of the SC coherent 
peak. It should be noted that the spectral weight in the 
antinodal region, observed by Kondo et al.) 3 - is consid- 
erably greater than that corresponding to simple cos 2^ 
dependence as shown schematically in Fig. 5a. 

The coherent peak disappears in the PG state above 
T c where spectral weight W PG ((f)) is zero in a broad angle 
range that can be referred to the nodal region (Fig. 5b). 
In the antinodal region, W PG (<p) increases rapidly up to 
a maximum when tf> — ► 0. Insulating order, which might 
be invoked to explain both W GP (<p) and W PG (<f), is not 
discovered for now. We believe that a competition of such 
a hidden order— with superconductivity is not the only 
qualitative explanation of observed spectral properties in 
the antinodal region. We have shown that K - pairing 
concept^ leads to a consistent explanation of the origin of 
the SC and PG states: K - pairing in the antinodal region 
gives rise to both these states. A decrease of W GP ((j)) at 
can be associated with non-trivial zero lines of 
the SC order parameter A sc (fc) shown in Fig. 3 that is 
can be explained in just the same way as a decrease of 
W CP (4>) at <fi — > 7r/4 due to d-wave node. It is evident 
that such angle dependence of the spectral weight of the 
coherent peak, appearing in both nodal and antinodal 
regions owing to different microscopic mechanisms of SC 
pairing, can be considered as a direct indication of the 
biordered SC state. On the contrary, since the PG state 
is associated with QSS wave function, one can expect 
that W PG (4>) — in the nodal region whereas nonzero 
W PG (4>) in the antinodal one is compared with W CP (<j>) 
that appears there below T c . Due to a random phase of 
the wave function of QSS and corresponding gap function 
A pg (k), nodal lines of these functions, that could be ap- 
parent in the antinodal region, cannot give a detectable 
contribution into a decrease of the spectral weight at 
(j) — ► 0. Expected angle dependencies of W CP ((f>) and 
W PG (4>) are shown schematically in Fig. 5. 

We believe the antinodal region with pronounced nest- 
ing of the FC gives rise to If -pairing whereas conven- 
tional pairing with K — develops in the nodal re- 
gion where the FC shows no signs of nesting. Thus, we 
do not oppose If-pairing with the conventional pairing: 
these two SC pairing channels with slightly overlapped 
domains of definition in the momentum space form indi- 
visible biordered SC state together. A passage from the 
antinodal region into the nodal one is accompanied with 
a redistribution of the spectral weight between these two 
pairing channels. 




<> tt/4 cf) <|> 71/4 <|> 



a b 

FIG. 5: Angle dependence of the spectral weight (schemati- 
cally), a. Spectral weight of the coherent SC peak (dashed 
lines 1 and 2 correspond to d-wave order parameter (~ cos 2(j>) 
and K -pairing order parameter, respectively, b. Antinodal 
spectral weight of the PG state. Here, <j>a is a conditional 
angle boundary separating the nodal and antinodal regions. 

IV. QUASIPARTICLE INTERFERENCE 

A rise of the coherence in the system of the antinodal 
K -pairs below T c should inevitably lead to interference 
effects inherent in the SC state. Bogoliubov QPI appears 
due to mixing of quasipatricle states with high spectral 
weight that results in a modulation of the local density 
of states (LDOS) in the real space. Such states, at given 
quasiparticle energy ((¥]), are disposed in vicinities of the 
points corresponding to maximal curvature of the isolinc 
E(k x ,k y ) = E = const. In the case of biordered SC state, 
pairing with zero total momentum dominates K -pairing 
in the nodal region, therefore, rj(k) = in this region 
due to the fact that s(—k) — e(k). Thus, quasiparti- 
cle spectrum in the nodal region turns out to be fully 
symmetrical with respect to the Fermi level. 

On the contrary, If -pairing dominates the pairing with 
zero momentum in the antinodal region including near 
rectilinear segments of the isolines in relatively small 
vicinity of the FC which primarily forms the singularity 
in the self-consistency equation. Because e(K/2 + k) « 
e(K /2 — k) in this vicinity, quasiparticle spectrum in the 
antinodal region should be slightly asymmetrical with re- 
spect to the Fermi level. Due to the fact that there is a 
considerable increase of the deviation from mirror nest- 
ing in a vicinity of conditional boundary separating nodal 
and antinodal regions, one can expect a pronounced in- 
crease of the asymmetry of quasiparticle spectrum in this 
vicinity. Such a statement is compatible with photoemis- 
sion data presented by Yang et al. 30 Also, it shows that 
an insulating state, competing with the SC one, cannot 
be considered as the only origin of the asymmetry ob- 
served by Yang et al»2£ Indeed, quasiparticle spectrum in 
the case of insulating state with gap function D(P) has 
the form 

E Q (p) = £(P) ± VV 2 (P)+D(P) 2 , (16) 

where 2f(p) = e(p) +e(p + Q), 2rj(p) = e(p) -e(p + Q). 
Therefore, imperfect nesting, that is a deviation from the 
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FIG. 6: Main scattering momenta determining nodal (a, in 
accordance with the "octet" model^) and antinodal (b) QPI 
pattern. FC and nodal line of the gap function are presented 
as solid and dashed lines, respectively. Small 'bananas' are 
showed as shadowed ovals. Antinodal sectors of the Brillouin 
zone are shadowed. Integers present the subscript enumerat- 
ing main scattering momenta. 

nesting condition (|12p . just as imperfect mirror nesting, 
results in the term (£(p) or r)(p) before the square root) 
that originates electron-hole asymmetry in both cases. 

To study antinodal quasiparticle spectrum qualita- 
tively, one can neglect rj(k) in Eq. (fH))) owing to the 
fact that antinodal segments of the FC appear as near 
rectilinear. For this reason, at low quasiparticle ener- 
gies, isoline shape can be analyzed in general form. In 
such a case, isolines enclose the gapless points of inter- 
section of the FC and the nodal line. These singular 
points can be found from equation system ^ K {k x , k y ) = 0, 

The nodal part of the FC can be approximated by an 
arc of a circle whereas the nodal lines of d - wave su- 
perconductor are straight lines k y = ±k x . As a result, 
the quasiparticle isolines become apparent as 'banana'- 
like closed curves^ Because Fermi energy e F exceeds 
d - wave SC gap magnitude A m considerably, e F ;» A, m , 
'banana' turns out to be very thin so that exactly its 
end points correspond to maximal curvature of the iso- 
line. This directly leads to the octet model of QPI in 
the nodal regional which defines a set of wave- vectors k t 
(1 < i < 8), corresponding to such end points, that de- 
termine LDOS pattern at given quasiparticle energy E. 
The octet model is presented in Fig. 5a where we define 
the main scattering momenta as q { = k i — k x . It should 
be noted that our definition of q i is somewhat different 
from that given by Kohsaka et el^ These two definitions 
are mutually complementary in the reciprocal lattice. 

As one can see from Fig. 3, singular points of antinodal 
quasiparticle spectrum E(k) are symmetrically, with re- 
spect to the antinodal directions, disposed on near nested 
pieces of the FC. To study a shape of the isolines, one 
can approximate these pieces by straight lines and any 
of the nodal lines in a small vicinity of the singular point 
by an arc of a circle. Evidently, under the condition that 
e F A m , all isolines in this vicinity are 'bananas' en- 
veloping a rectilinear part of the FC. Indeed, if E <C e F , 



an isoline of the the quasiparticle dispersion can be writ- 
ten in the form 

k x -k\ = -k\ ± a^Jn 2 E - n 2 y , (17) 

where n x — k x /k F , n y = k y /k F are dimensionless com- 
ponents of the relative motion momentum, k F is the 
Fermi momentum in the antinodal direction, n is di- 
mensionless radius of the nodal line, k e = E/2e F , 
a = 2e F /A m (k F a) 2 . As follows from Eq. ([17]), there 
are closed isolines only under condition that — k e < 
K y < k e- Therefore, a transversal, with respect to 
the FC, size of the isoline equals n t = 2k e <C 1. If 
E — > 0, closed isolines shrink into two singular points 
(±k ,0). A longitudinal size can be estimated as — 
a/kq + otK E — -\/ K o — aK E , therefore, closed isolines ap- 
pear in quasiparticle energy range < E < A m (k Q a) 2 . 
In addition, one can examine that K t -C k ; . Energy 
E m = A m (k Q a) 2 corresponds to a topological transition 
from closed, at E < E m , to opened, at E > E m , isolines. 
Because opened isoline has no points of considerable cur- 
vature, the topological transition should result in a degra- 
dation of the interference pattern. Contrariwise, due to 
a large curvature of the closed isoline in small vicinities 
of its end points, exactly these vicinities should primarily 
contribute into the QPI. Therefore, following McElroy et 
one can introduce a set of momenta q^E) connect- 
ing different end points. Here, subscript i runs from 1 to 
2n — 1 where n is a number of singular points of quasi- 
particle dispersion (j 16|) . Such main scattering momenta, 
defined as q i = k t — k x for any i ^ 1, should determine 
the real-space interference pattern. The pattern caused 
by the antinodal QPI turns out to be considerably more 
complicated in comparison with the nodal one even in 
the simplest case corresponding to the only closed nodal 
line in each of four crystal equivalent parts of the antin- 
odal region as shown in Fig. 6b. The full set of momenta 
q^ following from non-trivial momentum dependence of 
the gap function shown in Fig. 3, should result in the 
real-space antinodal QPI pattern that can be considered 
as originating from fairly uniform distribution of scatter- 
ing momenta. Therefore, it seems highly probable that 
the antinodal QPI pattern should be considerably more 
smooth with respect to the nodal one. 

It is clear that due to expansion of the closed isolines 
E(k) = E with an increase of E, there is a variation (ro- 
tation and decrease or increase of the absolute value) of 
the main scattering momenta, q i = q^E) with E. All 
of the nodal scattering momenta are dispersive, varying 
with E in accordance with the octet models On the con- 
trary, among the antinodal scattering momenta, there 
are some non-dispersive, such as "immobile" q 5 and q 14 
shown in Fig. 6b, that, owing to nesting feature of the 
FC, remain independent of E at small quasiparticle ener- 
gies. Strictly speaking, only such "immobile" q i can con- 
tribute into checkerboard real-space modulation in the 
SC stated 

Nodes of the antinodal quasiparticle spectrum result in 
the fact that, at a finite temperature, thermal equilibrium 
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quasiparticles are excited not only near d - wave nodes 
in the nodal region^ but in the antinodal one as well. 
Moreover, the equilibrium population of the antinodal 
quasiparticles may considerably exceed their population 
in vicinities of the d - wave nodes. This may occur if line 
of zeroes of the antinodal gap function, in contrast to the 
nodal one, turns out to be close to the FC in its extended 
vicinity. 

V. CONCLUSION 

Our concept of -Repairing in the cuprates is based on 
two complementary statements: 1° screened Coulomb re- 
pulsion is the underlying SC pairing interaction; 2° large 
momentum of JtT-pair arises due to nesting feature of the 
FC. SC ^-pairing, prevailing in the antinodal region of 
the momentum space, leads directly to uniform explana- 



tion of spatial inhomogeneity of both SC state in the form 
of checkerboard PDW and striped PG state formed by 
incoherent i^-pairs. AT-pairing, together with the con- 
ventional SC pairing with zero momentum prevailing in 
the nodal region, results in an indivisible biordered SC 
state which naturally explains the peculiarities of the an- 
gle dependence of the spectral weight both below and 
above T c . Complicated momentum dependence of the 
gap function in the antinodal region should lead to fairly 
reach antinodal QPI resulting in relatively smooth real- 
space interference pattern. 
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